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Abstract 

The aim of our paper is to formulate and solve problems concern¬ 
ing linear multiple periodic recurrence equations. Among other things, 
we discuss in detail the cases with periodic and multi-periodic coef¬ 
ficients, highlighting in particular the theorems of Floquet type. For 
this aim, we find specific forms for the fundamental matrix. Explic¬ 
itly monodromy matrix is given, and its eigenvalues (called Floquet 
multipliers) are shown. The Floquet point of view brings about an 
important simplification: the initial linear multiple recurrence system 
is reduced to another linear multiple recurrence system, with con¬ 
stant coefficients along partial directions. The results are applied to 
the discrete multitime Samuelson-Hicks models with constant, respec¬ 
tively multi-periodic, coefficients, in order to find bivariate sequences 
with economic meaning. 
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1 Discrete multitime multiple recurrence 

A multivariate recurrence relation is an equation that recursively defines a 
multivariate sequence, once one or more initial terms are given: each further 
term of the sequence is defined as a function of the preceding terms. Some 
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simply defined recurrence relations can have very complex (chaotic) behav¬ 
iors, and they are a part of the held of mathematics known as nonlinear 
analysis. We can use such recurrences including the Differential Transform 
Method to solve PDEs system with initial conditions. 

In this paper we shall continue the study of discrete multitime multiple 
recurrence, giving original results regarding Floquet theory of linear multiple 
periodic recurrences. As an example of application in economics, we for¬ 
mulate and solve a Samuelson-Hicks multitime multiplier-accelerator model. 
Such problems remain an area of active current research in our group. The 
scientihc sources used by us are: general recurrence theory n, 0,0, m, 
[21], [23], our results regarding the diagonal multitime recurrence [I], |5], Flo¬ 
quet theory 0, 0. i, CDl, 1221 and multitime dynamical systems jl2j-|2Uj. 

2 Linear multitime multiple recurrence with 
multi-periodic coefficients 

Let m > 1 be an integer number. We denote 1 = (1,1,...,!) G Z™. Also, 
for each a G {1,2,..., m}, we denote 1 q = (0,..., 0,1, 0,..., 0) G Z'", i.e., 
Iq, has 1 on the position a and 0 otherwise. 

On Z™', we dehne the relation “ < for t = ..., t™), s = (s^,..., s™), 

s <t if s" < t", Va G (1,2,..., m}. 

One observes that “ < ” is a partial order relation on Z™'. 

Let us formulate a theory similar to those of Floquet which brings im¬ 
portant simplihcation: in conditions of multiple periodicity or periodicity, a 
linear multitime multiple recurrence system 

x{t + Iq) = Aa{t)x{t), yt > to, Va G {1, 2,..., m} 

is reduced to a linear multitime multiple recurrence system 

y{t -F Iq) = B^y{t), Vt > to, Va G (1, 2,..., m}, 

with constant coefficients. 

We denote by Z one of the sets Z™- or {t G Z™-11 > ti} (with ti G Z™). 
Consider the functions Aq : Z —)■ Al„(iC), a G {1, 2,..., m}, which define 
the linear homogeneous recurrence 

x{t + la) = Aait)x{t), Va G (1,2,... ,m}, (1) 

with the unknown function x: {tG2^|t>to}—t iV"" = Ain,iiK), to G Z. 

We recall some notions and results of work [6] : Theorem [U Propositions 
[H [21 [3] and Dehnition [H 


2 


Theorem 1. a) If, for any (to, Xq) G ZxK^, there exists at least one function 
x: {tG2^|t>to}—)■ K^, which, for any t > to, verifies the recurrence dlD 
and the condition x(to) = xq, then 

Aa{t + l^)A^(t) = Ay{t + la)Aa{t), (2) 

'it e Z, Va,/3 e {1, 2,..., m}. 

b) If the relations ([2]), are satisfied, then, for any {to,xo) E Z x K^, there 
exists a unique function x: [t E Z \ t > to} ^ K"^, which, for any t > to 
verifies the recurrence ([T]) and the condition x{to) = xo- 

c) Let us suppose that Z = Z™' and that, for any a G {1,2, 

and any t E Z™, the matrix Aa(t) is invertible. If the relations ([2]), are 
satisfied, then, for any pair (to, Xq) E Z'" x K^, there exists a unique function 
x: Z'" —)■ K"^, which, for any t E Z™, verifies the relations ([I]), and also the 
condition x(to) = Xo- 

Proposition 1. Suppose that the relations o hold true. 

For each to E Z and Xq G fi4n(K) there exists a unique matrix solution 

X '. ItG^lt^toj" —t Xin(hl-) 

of the recurrence 


X(t + la) = Ao,(t)X(t), Va G {l,2,...,m}, (3) 

with the condition X(to) = Xq. 

For each to E Z, we denote 

x{' Ao)'- {t E Z \ t >to} ^ M.n(K), 

the unique matrix solution of the recurrence ([3]), which verifies X(to) = In- 

Definition 1. Suppose that the relations o hold true. The matrix function 

X(- ,■)■ {(t,s) E Z X Z\t> s} ^ Mn(K) 

is called transition (fundamental) matrix associated to the linear homoge¬ 
neous recurrence (fTj). 

For a E (1,2,..., m} and fc G N, we define the function 


Ca,k(t) 


Z ^ Mn(K), 

k 

Y\_Ao,(t + (k - j) ■ la) if k>l 

In if k = 0 . 


(4) 
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Proposition 2. Suppose that the relations ([2]) hold true. 

The matrix functions x(-) and Ca^ki.-) have the properties: 

a) x(t, 'S)x(s, r) = xif, ^), Vt, s,r E Z, with t > s > r. 

b) x{s,s) = In, Vs e Z. 

c) x{t + k ■ la, s) = Ca,k{t) ■ x{t, s), Vfc G N, Vt, s E Z, with t > s. 

d) Ca,k{t) = x{b + k- la, t), Vfc e N, Vt G Z. 

e) Ca,k{t + p- lp)Ci 3 ^p{t) = Cp^p{t + k ■ la)Ca,k{t), ^k.p G N, Vt G Z. 
/) For any t,s E Z with t > s, we have 


m—1 


X{t, s) = JJ C", 




a=l 


Y j.Q;+l j.ri 

,6 ? •••? 


■)-a 


,S 


^m—1 


g) For any t,s E Z, with t > s, the fundamental matrix xif, s) is invert¬ 
ible if and only if, for any a E {1,2,..., m} and any t E Z, the matrix Aa(t) 
is invertible. 

h) For any a E (1,2,..., m}, any k E N and for any t E Z, the matrix 
Ca,k{t) is invertible if and only if, for any a G {1, 2, ..., m} and any t E Z, 
the matrix Aa{t) is invertible. 

i) If'ia G (1,2, ...,m}, 'it E Z, the matrix Aa{t) is invertible, then 
it, s, to E Z, with t > s >tQ, we have xif, s) = xif, to)x{s, to)~^. 

j) Ifia G {1, 2, ..., m}, the matrix functions Aa{-) are constant, then 

C^,kit) = A^, ik EN,itE Z™, Vo G {1, 2,..., m}, 

x{t, s) = A^ ^ • ... • it, s E Z™, with t > s. 

Proposition 3. We consider the functions A ^: Z -E Ain{K), a E {1,... , m}, 
for which the relations ([2]) are satisfied. Let (tQ,xo) E Z x . Then, the 
unique function x: = A4n,iiK), which, for any t > to, 

verifies the recurrence and the initial condition x{to) = xq, is 


x{t) = x{t,to)xo, it>to. 


Ifia E (1, 2,..., m}, the matrix functions Aa{-) are constants, then 


x{t) = Af 


T-th) Ah^-to) 


..■Air-^o)x„ it>to. 


(5) 


2.1 Case of multi-periodic coefficients 

Definition 2. A function f: Z ^ M is called multi-periodic if there exists 
(Ti, T 2 ,..., Tm) E N™ \ {0} such that 

fit + Ta-la) = fit), itEZ, Va G (1,2,... ,m}, 

i.e., Ti ■ li, T 2 • I 2 , ..., Tm ■ Im are periods for the function f. 
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Proposition 4. Let T = (Ti, T 2 ,..., T^) G T 7 ^ 0. 

tKe consider the matrix functions Aa'. Z ^ AA^^K), a G {1,2,..., m}, 
for which the relations ([2]) are satisfied. Suppose that, Va, {3 G {1,2,..., m}, 
Vt G Z, we have Aa{t + Tp ■ 1^) = Aa{t). Then 

a) Ca,k{t + Ti 3 ■ 1^) = Ca,k{t), ya,j3 G {l,2,...,m}, Vt G Z; 

b) Ca^Tc{t)Ci3,Tp{t) = Ci3^Tpit)Ca,Tc,it), Vt G Z; 

c) x(t + Ta ■ la,s) = x(t,s) •C„,T„(s), Vt, s G Z, with t > s. 

Proof, a) It follows from the definition of Ca^k{')] this is either the constant 
fnnction or a prodnct of mnltiperiodic matrix fnnctions. 

b) According the step a), we have 

Ca,Ta{t + TjS ■ 1^) = Ca,Tc,{t), Cp^Tpit + Ta ■ Iq-) = Cp^Tp{t). 
According the Proposition [21 e), we have 

Ca,Tcii + Tp ■ li3)Ci3^Tp{t) = CjS^Tpit + Tq • 1q-)Cq,, Tc (t) • 

It follows that Ca,TA^)Cf}^Tp{t) = Cis^Tp{t)Ca,TSt)- 

c) Fix a. Fix s E Z. Let Fi, I 2 : {t G 11 > s} —)■ Ain{K), 

Yi{t) = x{t + Ta-la,s), Y2{t) = x{t,s) ■ Ca,Tc.{s), Wt > S. 

For each /3 G {1, 2,..., m}, we have 

Yiit+lp) = x(t + TQ-la + l^, s) = Ap{t + Ta-la)x{l'3-TaAa, s) = Ap{t)Yi{t)] 

Y2{t + lf}) = x(t + l/3,s)-C'a,T,(s) = A^(t)x(t,s) • C'„,T,(s) =A^{t)Y2{t). 

Conseqnently the fnnctions Fi(-) and F 2 (-) are both solntions of the recnr- 
rence ([3]). 

According the Proposition [H d), we have x(s + T„ ■ 1q,,s) = Ca,Ta,{^)j 
which is eqnivalent to x(s + ■ !«, s) = x(s, s)Ca,Ta{^)y i-®-) ^i('5) = ^ 2 ( 5 ). 

From the nniqneness property (Proposition [T]) it follows that Fi(-) and Y 2 {-) 
coincide; hence x(t + Ta ■ Iq, s) = x(t, s) ■ Ca,Ta{s), Wt > s. □ 

It proves withont difficulty the following result. 

Lemma 1. Let A, B E Al 2 (C), such that A 7 ^ A/ 2 , VA G C. Then 

AB = BA if and only if there exist z,w E C, such that B = ZI 2 + wA. 
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Proposition 5 . Suppose that the matrices Pi,P2,...,Pm G A^„(C) com¬ 
mute, i.e., PaP/s = PjsPa, Vo, /9 G {1, 2, ..., m}. 

Let us assume that one of the statements i) or ii) below is true. 

i) For each a G {1,2,..., m}, the matrix Pa is diagonalizable. 

ii) n = 2 and for any a G (1, 2,. .., m}, the matrix Pa is invertible. 

Let ki, k2, ... ,km G M*. Then, for any a G {1,2,..., m}, there exists 
Qa £ A 1 „(C) such that 

a) = Pa, Vo G {1, 2, ..., m], 

b) QaQy = QyQa, Va,/3 G { 1 , 2 ,.. . ,m}. 

Proof. For di, d2, ■ ■. ,dn G C, we denote diag((ii, d2, ■ ■ ■, dn), the matrix 
from Mn(C), with main diagonal elements di, d2, ... ,dn (in this order) and 
in rest 0 . 

If Hypothesis i) is true and if PaPy = PyPa, then the matrices Pa 

are simultaneously diagonalizable, i.e. there exist T G Aln(C), T invertible 
and there exist Xa,2, ■ ■ ■ , Xa,n^ C such that 


Pa = T ■ diag(Ao,i, Aa, 2 , • • • , Aa,n) ' T \ Vo G {1,2,..., m}. 

There exist 9ay G C, such that = A^j, Va G {1,2, ...,m}, Vj G 
{l, 2 ,...,n}. 

For Qa = T ■ diag( 6 'Q,^i, 2 ,..., Oa^n) ■ T~^, a G {1,2,..., m}, easily hnds 
that statements a) and b) hold. 

Suppose that the hypothesis ii) is true. 

If Vo, 3Aa G C, such that Pa = Xah, then the matrices Pa are diago¬ 
nalizable and hence the hypothesis i) is true, and this case has already been 
treated. 

We left to study the case in which at least one of the matrices Pa is not 
of the form Pa = Aa/ 2 , with A^ G C. After a possible renumbering we can 
assume that this matrix is Pi. Hence Pi 7 ^ A/ 2 , VA G C. 

Since Pi Pa = PaPi, from the Lemma[I]follows that there exist Za, Wa G C, 
such that Pa = Zal 2 + 'UJaPi (Vo G {1,2,..., m}). 

If Pi is diagonalizable, then the matrix Zah + WaPi is also diagonalizable. 
So assuming i) is satished, and this case has already been treated. 

If the matrix Pi is not diagonalizable, then there exists a matrix T G 


M.2{P.), T invertible and there exists A G C, such that Pi = T ^ ^ ^ J 
Consequently, we have Pi = XI2 + S, where S' = T f ^ ^ | T~^ 




We observe that Pq = (za + Xwa)l2 + WaS, Va G { 1 , 2 ,..., m}. Since 
S^ = O2, we deduce that the matrix S is not invertible; hence Za + Xwa 7^ 0 . 
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It follows that there exists G C, 7 ^ 0, such that = Za + Xwa- Let 

_ Wq, 

kaU^°‘~^ ' 

We choose the matrix = Uah + VaS. Since = O 2 , Vj > 2, we 
obtain 

+ kaU’^°‘~^VaS = {Za + XWa)l2 + WaS = Pa, Vtt G {1, . . . , m}. 
The equality QaQis = Q/sQa is obvious. □ 

Proposition 6 . Let to G IP, fixed. We denote Z := \^t E |t>to}- 
LetT=(Ti,T 2 ,...,T^)GN"*, T^O. 

Consider the matrix functions Aa - Z ^ W1„(C), a G {1, 2,, m}, for 
which the relations o are satisfied. Suppose that, \/a,(3 G {1,2, 

Vt G Z, Aa{t + • Ijf) = Aa{t) and Aa{t) is invertible. 

We denote <h(t) := x{tAo), t E Z, and Cq = Ca,Ta{to)- 

Let Fi = {a G {1,... ,m} | > l}, F 2 = {a G (1,... ,m} | = O}. 

For each a E Fi, we choose Ba E Aln(C) such that Bff^ = Ca- 
For each a E F 2 , we choose Ba = In- 

If for any a, (3 E Fi, we have BaBp = BpBa, then there exists a function, 
P\ Z ^ M-niC), such that 

P{t + Ta ■ la) = P{t), Vt E Z, Va G {1, 2,..., m}, and 
^t) = P{t)BtBf-...-BZ, Vt>to. 

Proof. We remark that for any a, (3 E (1,2,..., m}, we have BaBp = BpBa- 
If a G F 2 , i.e. Ta = 0, then Ca = In- We observe that the equality 
5^“ = Ca is true. Hence, for any a G (1, 2,..., m}, we have Cq = 5^“. 

Let a E Fi, i.e. Tq > 1. Since the matrix Cq is invertible (Proposition [2]), 
from Bff^ = Cq it follows that the matrix Ba is invertible. 

Hence, for any a G (1, 2,..., m}, the matrix Ba is invertible. 

We dehne the function P\ Z ^ Aln(C), 

P{t) = ^{t)BfBf -...- Bff"', Wt > to. ( 6 ) 

Since from the formula ([ 6 ]) it follows immediately the equality 

Ht)=P{t)BfBf-...-BZ, 

it is sufficient to show also that P(-) is a multiperiodic function: 

P{t + Ta • 1q) = x{t + Ta • 1q, to)BfBf -...- • ... ■ Bff^. 


7 



According the Proposition HJ c), we have 


Xit + Ta ■ la, to) = x{t,to) ■Ca,TAto) = ^{t)Ca = 

We obtain 

Pit + Ta ■ 1„) = mBl-Bfpf ■ ... • • ... ■ 5-*™ = 

= mBfBf • ... ■ 5-*“ ■ ... • 5-*”* = Pit). 

□ 


Theorem 2. Let to G Z™, fixed. We denote Z := \^t E Z™' 
LetT=(Ti,T2,...,T^)eN-, T^O. 

Consider the matrix functions Aq : Z —)■ Ad„(C), a G {1, 2,, m}, for 
which the relations ([2]) are satisfied. Suppose that, \/a,/3 G {1,2, 

Vt G Z, Aait + • Ifi) = Aait) and Aait) is invertible. 

We denote $(t) := xit,to), t E Z. 

Let us assume that one of the statements i) and ii) below is true, 
i) n = 2. 

n) Va G (1, 2,..., m}, Wt G Z, the matrix A^t) is Hermitian and 

Aait)Aait + k-la) = Aait + k ■ la)Aa(t), Vt G Z, V/c G N, Vo G {1, . . . , 777.}. 

Then there exists a function, P: Z ^ A1„(C), and also there exist the 
constant invertible matrices Bi, B 2 , .. ., Bm E Aln(C), such that 
a) Pit + Ta ■ la) = Pit), yt E Z, Va G {1, 2,..., 777 }; 
h) BaBfi = BfiBa, Va,/?G{1,2,...,777}; 

c)m = P{t)BtBt -...-bz, yt>to. 

Proof. Let Ca = Ca^rdto)- The matrices Ca are invertible (Proposition [2]). 

One observes that in the hypothesis ii), the matrices Ca are Hermitian, 
hence diagonalizable. 

Let Fi = {a G (1,..., 777 } | > l}, F 2 = {a G {1,..., 777 } | Tq, = O}. 

For each a E F 2 , i.e. Ta = 0, we choose Ba = In- 

For the set of invertible matrices {Oa | a G Fi}, we apply the Proposi¬ 
tion [5l The Proposition m says that CaCp = CpCa- Note that either n = 2 
or the matrices Ca are diagonalizable. Hence the hypotheses of the Proposi¬ 
tion [5] are trne. Hence, for each a G Fi, there exists Ba E A1„(C) snch that 
Bffi = Ca and 

BaBp = BpBa, Va,/1 G Fi. 

From Bffi = Ca and Ta > 1, it follows that the matrix Ta is invertible. 


Since for a G F 2 , we have Ba = Im h follows 


= Va,/? e 

We observe that the hypotheses of the Proposition [B] are true. Conse¬ 
quently, it follows automatically the points a) and c). □ 

Theorem 3. Let Iq G Z™, fixed. We denote ^ := G Z"* 
LefT=(Ti,T2,...,TjGM-, T^O. 

Consider the matrix functions Aa'. Z ^ W1„(C), a G {1, 2,, m}, for 
which the relations ([2]) are satisfied. Suppose that, Va G {1,2,..., m}, Aa(t) 
is invertible, \/t G Z. We denote <h(f) := xifAo), t & Z. 

We assume that there exists a function P: Z ^ A1„(C) and there exist 
the constant invertible matrices Bi, B 2 , ..., B^ G A1„(C), such that the 
relations from the steps a), b), c) of Theorem^ to be satisfied. 

We consider the recurrences 

x{t + la) = Aait)x{t), Vf > to, Va G {1,2,... ,m}; (7) 

y{t + la) = BaVit), Vt > to, Vtt G {1,2,... ,m}. (8) 

If y(t) 'i'S a solution of the recurrence ([8]), then x(t) := P{t)y(t) is a 
solution of the recurrence o. And conversely, if x(t) is a solution of the 
recurrence dZD, then y(t) := P(t) ^x(t) is a solution of the recurrence ([H]). 

Proof. Since the matrices Ba commute, it follows that the recurrence ([8]) 
has the existence and uniqueness property of solutions (see Theorem [1]). 

The matrix <h(t) is invertible (Proposition [2]). From the equality of the 
point b) in the Theorem [21 it follows that the matrix P{t) is invertible. 

Let y{t) be a solution of the recurrence ([8]) and x{t) := P{t)y{t)] hence 
y{t) := P{t)-^x{t). 

y{t + la) = Bay{f) P{t + la)~^x{t + !«) = BaP{t)~^x{t) 
x{t + lo) = P{t + la)BaP{t)~^x{t). 

^ x{t + la) = $(f + la)Bf^'Bf^^ ■ . . . • • . . . ■ Bfi!^ ■ BaP{t)-^x{t) 

^ X{t + la) = Aammf^'Bf^^ ■ . . . • 5 -*“ • . . . ■ Bfi!^ • P{t)-Ht) 
x(t + la) = Aa(t)P(t) ■ P(t)~^x(t) x(t + la) = Aa(t)x(t). 

Like it proves the converse. □ 
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Conjecture 1. Suppose that the invertible matrices Pi, P 2 ,..., Pm £ 
commute, i.e., PaP/s = P/sPa, Va, /9 G {1, 2,..., m}. 

Let ki,k 2 , ■ ■ ■ ,km G N*. Then, for any a G there exists 

Qa G A4n(C) such that 

a) = Pa, Va G {1, 2,..., m}, 

b) QaQy = QyQa, Va,/3 G {1,2 ,... ,m}. 

Conjecture 2. Let to G TP, fixed. We denote := {t G TP 1 1 > to}- 
LetT={Ti,T2,...,Tm)eN^,T^0. 

Consider the matrix functions Aa - Z ^ Aln(C), a G {1, 2,..., m}, for 
which the relations ([2]) are satisfied. Suppose that, \/a,/3 G {1,2, ...,m}, 
Vt G Z, Aa{t + Tp ■ \jfi) = Aa{t) and Aa{t) is invertible. 

We denote <l>(t) := xifAo), t E Z. 

Then there exists a function, P: Z ^ Ain{C), and also there exist the 
constant invertible matrices Bi, B 2 , ..., Bm G A^„(C), such that 

a) P{t + Ta ■ la) = P{t), 'it E Z, Va G {1, 2,..., m}; 

b) BaBy = ByBa, ia, fi E {1,2,... ,m}; 

c) m = pmtBt-...-Bf:, vt>to. 

Remark 1. From the Conjecture 1 it follows the Conjecture 2. 

The proof is the same as that of Theorem [2l One uses the Conjecture 1 
instead of the Proposition [5] to show that the hypotheses of the Proposition [6] 
are satisfied. Then one applies Proposition^ 


2.2 Example and commentaries 

Example 1. Let us consider m = 2, n > 2, T 2 E T 2 > 2, 


Q = 



cos ■ 


sin 


TT 

71 

To. 


. VT \ 
— Sin — 
T 2 

71 

COS — 

T 2 


Si = 


1 0 
0 -1 ^ ’ 




02,n-2 

In-2 

If n = 2, we set Q = Qi and S = Si. 


Si 

On- 2,2 


02,n-2 

In-2 


dacd n > 3. 


We have = Qj, Q-^ = , S^ = Si = Sj, S-^ = S = S^, 

and hence the matrices Qi, Q are orthogonal and Si, S are symmetric and 
orthogonal. 


Powers: Q? = -I 2 , QT 


I 2 , = 


-I 2 

\On-2,2 


02,n-2\ ^2T2 

In-2 


In. 
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Also, we have SiQi = and SQ = Q ^S. Hence Q ^ = SQS and 

g-^ = SQ’^S or SQ-'^ = Q'^S, Wk e Z. 

O'TT Oyr 

If SQ = QS, then Q~^ = = In cos — = 1 and sin — = 0, 

T 2 T 2 

what can not becanse T 2 > 2 . Hence SQ 7 ^ QS. 

Let us consider the matrices 


Ai,A2:Z2^M„(C), 


Ai{t\e) = Q^" SQ-^" = Q^^" S, A2{t\e) = Q, 

Since the matrix S is symmetric and the matrix Q is orthogonal it follows 
that the matrix is hermitian (it’s really real and symmetric). The 

matrix is also orthogonal. 

Unfortunately the matrix is not hermitian (symmetric). This is 

the only hypothesis of Theorem [2] which does not take. 

Let us compute Ai\ 



QfS^ 02,n-2\ ^2t^ 

On-2,2 In-2 ) ’ 


/ 27rt^ 

cos —z— 
T 2 

2Tlt^ 
sin — 7 — 


27rt^\ 

— sin —z— 

T 2 . 

2nt^ ’ 
cos —z— 

T 2 / 


, with Qfs, = 

V^n- 2,2 -'n -2 / 


/ 27rt^ 

cos —Z— 
T 2 

27rt^ 


sm 


27rt^ \ 
27rt^ 


sm ■ 


To 


cos ■ 


T 2 / 


A2{t^ + l,t^)Ai{t\t^) = QQ^'SQ-^" = Q^"+^SQ-^\ 

Hence Ai{A,t‘^ + 1 )^ 2 (^^ 5 ^^) = A 2 {A + 1, t^)Ai(t^, t^), i.e., the relation ([2]) is 


true. 

Obviously that Ai{A + fc,t^) = Ai{A + k,t^)Ai{A E N, 
V(t\ t^) G Z^, since Ai(-, •) is constant with respect to the hrst argument. We 
have also A 2 {A,t‘^)A 2 {A, + k) = A 2 {A, + k)A 2 {A, t^), Vfc E N, V(t^, t^) E 

I?. 


Now, obviously that AiiiS + l,f^) = Ai(U,f^), V(t^,t^) G 1?. AiiAA"^ + 
T 2 ) = = Q'^^^S, since g^'^^ = In- Hence Ai(t^,t^ + T 2 ) = Ai(t^,t^), 

V(U,t2) e ZT 

Since A 2 (-, •) is a constant function, it follows that 


A2{t^ + lA") = A2{t\t^),A2{t\t^+T2) 


A2{t\t^), eZ\ 
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Was also observed that Ai{t^,t‘^)A2{t^,t‘^) ^ A2(t^, t^), V(t^,t^), 

because if we had equality would result Q* SQ~^ Q = QQ^ SQ~* , equivalent 
to SQ = QS, what is false (we proved above that SQ ^ QS). 

According Theorem [H for any to = (tl,to) ^ and any xq G C^, there 
exists a unique solution a:: —)■ C" of the double recurrence 


x{A + l,f2) 

x{t^,t^ + 1) 


AiiA, t‘^)x{A,t‘^) 
A2{t\t^)x{t\t^), 




(9) 


which verifies the initial condition x(to) = xq- Let us show that this solution 
can be written in the form 


x{t\i^) = g*'^*'-*og-*oa:o, v(f\tg e zt 


( 10 ) 


Indeed, 

x{A + l,f) = 

A,{t\t^)x{t\t^) = = g*'5i+*'-*og-'oa:o. 

x{t\t^ + 1) = = A 2 {t\t^)x{t\f) 

and the equality x{tl,to) — is obvious. 

It follows that the fundamental matrix is x(t,to) = “*og-*o. We 

select to = ( 0 , 0 ); let <h(f^,f^) := x ((^^5 ^^); ( 0 , 0 )) = . 

We shall determine the matrices Bi, B 2 and P(-) as in Proposition [HI 
with Ti = 1. 

Cl = Ci^Ti{to) = C'i,i(0, 0) = Ai(0, 0) = S, 

— I 2 02.n- 


C2 = C2,TM(t)=Q^- = 


_ nT2 _ / “-'2 <-^2,^-2 

On-2,2 In-2 


We look for the matrices Bi,B2 G A 1 „(C) that satisfy B^^ = Ci, B^^ = 


C2, i.e., Bi = S, and B^^ = [ ^ 

\ C '„-2 9 Jr 


'n-2,2 J-n-2 

number z E C such that = _x^ select B2 = 


. Since there exists a complex 
ZI 2 02,n-2 


On- 


n—2,2 In—2 


and 


conclude that B^^ = 


On- 2,2 


02,n-2 

In-2 


We have also B1B2 = B2B1 (i.e., SB2 = B2S), since Bi and B2 are 
diagonal matrices. We see that the assumptions of Proposition [6] are satisfied 
(with Ti = 1). 

The relation <h(f\f^) = P(A,t^)Bi Hf i^ equivalent to 


gt-^n _ p(^t\C)BlBl^ ^ = P{t\f)Bl 


P{t\C) = Q^^B^^^ 
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_ ( QX 02,n-2 

On-2,2 In-2 


X-* T 

J -19 


On- 


n-2,2 


'n-2 

+2 „_+2 


02,n-2\ _ [ ^ Ql 02,n- 


On- 


n—2,2 


P(t\t" + T2) = 


On-2,2 


In- 


n—2 


= Pit\0) 


‘■n-2 


P((0 + Ti, = p(tl + = Qi p-* = 

f ,y-t^-T2(^t^+T2 f) \ 

p{t\e + T2) = (^ „ ). 

V On-2,2 ^n-2 J 

Since = —1 and Qp = —I 2 , it follows 

02,n-2\_ 


We verified that P{t^ + Ti,t^) = P(f^,f^), P(f^,+ T 2 ) = P(f^, P). 

Let’s check also the conclusion of Theorem [3l 

According the Proposition [3l any solution of the recurrence ([8]) is of the 
form = B^B^v, with v G C"". On the other hand, we have 

P{t\t^)y{t\t^) = Q^"b^^"bXbXv = Q^'bXv 


and we observe that x{P,t‘^) := P{P,t‘^)y{P,t‘^) is indeed a solution of the 
recurrence dH]), which verihes the condition a;(0,0) = v, according flTU]) . 


2.3 Case of periodic coefficients 

Definition 3. Let T G N™, T 7 ^ 0. The function f: Z ^ M is called 
periodic of period T if f {t + T) = f (t), yt e Z. 

Proposition 7. Let T = (T\T^,... G M™', T 7 ^ 0. We consider the 
matrix functions A^: Z ^ Ain{K), « G {1,2,... ,m}, for which the rela¬ 
tions ([2D are satisfied. Suppose that, Vo G ( 1 , 2 ,..., m}, the matrix function 
Aa{-) is periodic of period T. Let C(t) = x(t + T,t). Then 
0 ) Ca,k{.-) is periodic of period T; 
h) x{t + T,s) = x(t, s) ■ C{s), Vt, s E Z, with t > s; 
c) 0(s) = C'l,Tl(s^ s^+TX s"^+T^)C 2 ,t^{s\ sX s^+TX 

(s^, . . . , S™ s’" + T™’)C'm, S^, . . . , s”* ysEZ. 

Proof, a) follows from the dehnition of Ca,k{')''i this is either the constant 
function In, or a product of matrix functions, periodic with the period T. 
b) Fix s E Z. Let Yi,Y 2 : {t E Z\t > s] ^ M.n{K), 


Yiif) = Xif + T,s), Y 2 {t) = xit,s) ■ C{s), \ft>s. 
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For each a G {1, 2,, m}, we have 
Yi{t + Ic) = xit + T + la, s) = Aa{t + T)xit + T,s) = Aa{t)Yi{t); 

Y 2 {t + 1 „) = xit + la, s) ■ C(s) = Aa(t)x(t, s) ■ C(s) = Aa(t)Y 2 (t). 

It follows that the functions Fi(-) and Y 2 (-) are both solutions of the recur¬ 
rence dSj); we have also x(s+T, s) = x{s, s)C'(s), i.e., Fi(s) = ^ 2 ( 5 ). From the 
uniqueness property (Proposition [T]) it follows that Fi(-) and F 2 (-) coincide; 
hence x{i + T,s) = x(t, s) ■ C{s), Wt > s. 

c) follows directly from Proposition [21 /). □ 

Theorem 4. Let to G fixed. We denote Z := {f G Z™' 1 1 > to}- 
Let T = T™) G M"*, T ^ 0. 

Consider the matrix functions Aa - Z ^ j\4n(C), a G {1,2,... ,m}, for 
which the relations ([2]) are satisfied. Suppose that, Va G {l,2,...,m}, the 
function Aa{-) is periodica de perioadd T and Aa(t) is invertible, Vt G Z. 

We denote <I>(t) := x(t,to), t E Z. 

Then there exists a matrix function P: Z ^ j\4n(C), periodic of period 
T, and also there exist the constant invertible matrix B G Mn(C), such that 
<h(t) = Vt > to, where |t| = t™'. 

Proof. The matrices Ca,k{t) are invertible (Proposition [2]). From the Propo¬ 
sition [71 c), it follows that the matrices (^(s) are invertible. Hence there 
exists an invertible matrix B G A1„(C), such that = C'(to), where 

_j_ _j_ _j_ 

We dehne the function 

P-.Z^MniC), P(t) = Vt> to¬ 

ft is sufficient to show that P(-) is periodic of period T. 

P(t + T) =x(t + T,to)P-l'+'^l. 

According the Proposition [3, h), we have 

X(t -F T, to) = x(t, to)C'(to) = <F(t)pl^L 


We obtain 


P(t + T) = $(t)pl^lp-l*+^l = d)(t)p-l*l = P(t). 


□ 
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Theorem 5. Let to G fixed. We denote Z := G Z"^ 1 1 > to}- 
Let T = {T\ T2, ..., T^) G AA™, T ^ 0. 

Consider the matrix functions Aa'. Z ^ A4n(C), a G {1,2,... ,m}, for 
which the relations ([2]) are satisfied. Suppose that, Va G (1,2,..., m}, Aa(t) 
is invertible, Vt G Z. We denote ^(t) := xifAo), t E Z. 

We assume that there exists a matrix function P: Z ^ periodic 

of period T, and also there exist a constant invertible matrix B G A4n(C), 
such that <l>(t) = P{t)B^^^, Vt > to- 
We consider the recurrences 

x{t + la) = Aait)x{t), Vt > to, Va G {1,2,... ,m}; (11) 

y{t + la) = By{t), Vt > to, Va G {1,2,... ,m}. (12) 

U yif) 'Is a, solution of the recurrence ffT^ . then x{t) := P{t)y{t) is a 
solution of the recurrence cm- And conversely, if x{t) is a solution of the 
recurrence ca, then y{t) := P{t) ^x{t) is a solution of the recurrence ffT^ . 

Proof. The matrix ^{t) is invertible (Proposition |2]). 

From <F(t) = P(t)B^^^ it follows that the matrix P(t) is invertible. 

Let y(t) be a solution of the recurrence flT^ and x(t) := P{t)y{t)] hence 
y{t) := P{t)-^x{t). 

y{t + la) = By{t) P(t + la)~^x{t + la) = BP{t)~^x{t) 

x{t + la) = P(t + la)BP{t)~^x{t). 

^ x{t + la) = <h(t + la)P"^'*'+^^PP(t)"^a;(t) 
x{t + la) = Aa{t)^{t)B~^^^^^P{t)~^x{t) 

X{t + la) = Aait)P{t) ■ P(t)~^x{t) X{t + la) = Aa{f)x{t). 

Like it proves the converse. □ 

3 Discrete multitime Samuelson-Hicks model 

We assume that t = (t^,..., t™) G N™ is a discrete multitime. Having in mind 
the discrete single-time Samuelson-Hicks model m. we introduce a discrete 
multitime Samuelson-Hicks like model based on the following economical el¬ 
ements: (i) two parameters, the hrst 7 , called the marginal propensity to 
consume, subject to 0 < 7 < 1 , and the second a as decelerator if 0 < a < 1 , 
deeper if a = 1 or accelerator'll a > 1 ; (ii) the multiple sequence Y{t) means 
the national income and is the main endogenous variable, the multiple se¬ 
quence C{t) is the consumption] (iii) we assume that multiple sequences Y(t), 
C{t) are non-negative. 
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3.1 Constant coefficients Samuelson-Hicks model 

We propose a first order discrete multitime constant coefficients Samuelson- 
Hicks model as first order multiple recurrence system 


(yit + w\ 

\C{t + Ijs)] 


+ a 

V 7 



Vt e M”", V/3 e 


F(0) = Fo, C(0) = Co; 


with a, 7 G C \ {0}. 

f 'y y 

The matrix of this double recurrence is A = 

V 7 

the Example Proposition [31 we have the solution 


l'’' j. According 


(Yit) 

[C{t) 



(13) 


Let ri,r 2 be the roots of the characteristic polynomial — (7 + a)r + a of 
the matrix A. It proves easily by induction that whether ri 7 ^ r 2 , then 




r 2 — r\ ^ r 2 r\ — rir^ 
r2 - ri r2 - ri 


■I2 


Vfc e N, 


and whether ri = r 2 , then 


= kr’l-^A -{k- l)r^/ 2 , Vfc e N. 


(14) 


(15) 


We set k = |t| = + C + ... +in the formula flTT)) or fiT5|) and we 

obtain the relation flT3|l which gives concrete expressions for Y{t), C{t) in 
function of ri, r 2 . 


3.2 Multi-periodic coefficients Samuelson-Hicks model 

Let us use the variable parameters 

m 

7; [j{teZ^\t>-lp] 

P=i 

such that 7 (t) 7 ^ 0 and 7 (t) + a{t) ^ {0; 1}, Vf G N™. They define a discrete 
multitime multiple Samuelson-Hicks model 

(yc[t + 1 ^)) ^ ^ ^ 2 ,..., m}, 
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F(0) = Fo, C^(0) = Co. 

The matrix 

/ 7(t) + a(t) — 7 r^^\ 

1 VteN”^, v/?e{l,2,...,m}, (16) 
V 7(C 0 J 

must satisfy the relations ([2]), i.e. 


We denote x{t) = 


Aj^it + — A^{t + li3)Ap{t), 

'it G e {l,2,...,m}. 

y{t)\ 
cit) ■ 


( 17 ) 


The Samuelson-Hicks recurrence writes 


x{t + l/s) = Aj 3 {t)x{t), it e N™, V/3 e {1, 2,..., m}, 

a;(0) = xq. 

The relation f[T7)l is equivalent to 

{lit) + a(t))[j(t + If,) + a(t + If,)) - * 

{'j{t) + a{t))j{t + If,) * 

(7(t)+ a(t))(7(t + 1;3)+ a(t + 1^)) - * 

( 7 (t) + a(t)) 7 (t + I/ 3 ) 

It follows that 

7(t + 1^) = 7(f + I/ 3 ) 

and 

= ( 7 (i) + a{t))a(t + V) - 


( 18 ) 


( 19 ) 


By induction, one obtains 'j{t + k ■ If,) = 'y{t + A; • l_a), V/c G N, Vt G N™, 
V/x, V/3 and 

7(A) = 7((A\ A^ ..., 0) + • 1^) = 7((A\ A^ ..., A'"“\ 0) + ■ li) 

= 7(a' + A^ ..., 0) = 7((A' + r, ..., 0, 0) + ■ l™_i) 

= 7((A^ TA'", A^ ..., 0, 0) +r-i ■ li) = 7(f^..., 0, 0) 
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= ... = 7(|i|,0,...,0). 

where |t| means the sum + ... + t™'. 

Let f{k) := 7 (fc, 0,..., 0), fc G M U { —1}; we have obtained ^{t) = f{\t\). 
The relation flT^ becomes 

( 7 (t) + a{t))a{t + 1^) - a{t + 1^,) = {'y{t) + a{t))a{t + l/j) - a{t + 1^) 

if and only if 7 (f) + a{t) 7 ^ 1. Then a{t + 1^) = a{t + 1^); analogously to 
7 (-), it is shown that there exists g{k), k eN, such that a{t) = g{\t\). 

(f(]A q(]A \ 

We denote A{k) = ( f{k — 1) , /c G N. On the other 

V fik) 0 J 

hand, we have A/sit) = A{\t\), Vf G N™, V/5 G {1, 2,..., m}; and immediately 
notice now that, in this situation, the relationships ffT7|l are satished. 

We showed that the relations drzn are satisfied if and only if there exist 
the functions f : NU { — 1} —)■ C, ( 7 : N —)■ C, such that 

m 

7(«) = /(|i|), V«e a(«)=g(|«|), Vi £ N”. 

13=1 


We consider the recurrence 

z{k + l) = A{k)z{k), VfcGM. (20) 

Let 2 ;: N —)■ be the solution of the recurrence fl2UD . with 2;(0) = Xq. 

Let x: N™' —)■ C^, x{t) = z{\t\), \/t G M™'. Then 

x{t + 1 ^ 3 ) = z{\t\ + 1) = A{\t\)z{\t\) = Ag{t)x{t). 

Obviously a:(0, 0,..., 0) = z(0) = xq. Hence, the following result is true. 

If z: N ^ is the solution of the recurrence (120]), with 2:(0) = Xq, then 
the function 

x:N^ x{t) = z{\t\), VtGN™ 

is the solution of the recurrence flTSD . which verifies a;(0) = Xq. 

Let p G N and 

Op: Wl2(C), 

f TT A(k + p — i), if p > 1 
Cpik) = {l\ y p k P- 

1 /2, ifp = 0. 
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The solution of the recurrence (EOD, with 2 ;( 0 ) = Xq is z{k) = Ck{0)xo, 
V/c G N. Hence, the solution of the recurrence (ITS]) , which verihes a;( 0 ) = Xq, 

is 

x{t) = Cit\{0)xo, 

Let x('! •) the fundamental matrix associated to the recurrence flTSl) and 
let Xi(‘ 5 ‘) the fundamental matrix associated to the recurrence fl2U]) . We 
denote <h(f) = x(t,0), t G and \h(fc) = Xi{k,0), /c G N. 

From the above observations immediately yield \h(/c) = Ca;(0), V/c G N, 
and <F(f) = C|i|(0), Vf G N™'. Hence <F(f) = th(|t|), Vt G M™'. 

The foregoing first order homogeneous multiple recurrence is called multi- 
periodic if its coefficients are multi-periodic or if the matrix A^{t) is multi- 
periodic. A prominent role in the analysis of a multi-periodic recurrence is 
played by so-called Floquet multipliers. 

Let us consider a discrete multitime multi-periodic coefficients Samuelson- 
Hicks model. For that, suppose 3/r, 35 and 3T G N* such that = 

Afj_{t), \/t G equivalent to A{F -|-f^ . -t-f™ -|-T) = A{A -\-F .. . +f™'), 

Vf G N’". Set F = k and F = 0, for jd >2. One obtains A{k -\-T) = A{k), 
V/c G N. And here follows immediately that A^(f-|-T ■ Ip) = A(s{t), Vf G N’", 
Vp, V/3. 

Consequently, 

z/3p, 35 and 3T G M* a.i Afj_{t -|- T ■ I 5 ) = Ap(f), Vt G M™', then 
Ap{t + T • Ip) = Affit), Wt G Vp, V/3 
and the function A(-) is in fact periodic with the period T; which is equivalent 
to the fact that the functions f and g are periodic, of multi-period T. 

Let us suppose that the matrices A{k) are invertible; it appears the con¬ 
dition g{k) 7 ^ 0, V/c G N. 

Let us compute the matrices for fhe recurrences flTSD . 

fl2U]) . with to = 0. Obviously, in the case of the recurrence (I2UD . we have 
Tfs = T, V/3. 

For fITSD . we have a single such matrix, namely C = Ct{0) = T(0). 

For the recurrence (120]) . we have Cp = C^,r(0) = C't(O) = C. 

If the matrices A{k) are invertible, then C is invertible; hence, there exists 
B G A42{C), such that = C. Obviously, this is equivalent to BF = 
(hence = B, V/9). 

Hence, we can apply the Proposition [ 6 ] for the recurrence flTSl) (is in fact 
the case m = 1 ) and fl 20 |l (the matrices Bp commute, since they are equal). 
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Consequently, there exists R{k) and P{t) such that 


R{k + T) = R{k), VfceN, 

Vfc e N, 

P(t + T-1^) =P(t), VteN™, V/?, 
$(t) = VteN™. 


But <h(t) = ^(|t|) ^ P(t)pl*l = P(|t|)pl*l ^ P(t) = P(|t|). Hence, 
we have obtained: P{t) = P(|t|) and <h(t) = P(|t|)Pl*l, Vt G N™. 

Suppose that we are in case of multi-periodic recurrence of the type 
fllSp (discrete multitime multi-periodic coefficients Samuelson-Hicks), i.e., 
the functions / are g periodic, with the period T > 1 (equivalent to y4(-) 
is periodic with the period T). 

The matrix 

T 

C = Cg = Cf,,T,{0) = C'^,r(0) = C't(O) = 1[A{T - j) 

j=i 

is called monodromy matrix associated to the (multi-periodic) recurrence 

m- 

According the Proposition 01 c), we have T(A;-|-T-lyj) = T(fc)-C, Vfc G N. 
By induction, it follows T(fc -|-pT ■ I/ 3 ) = T(A;) ■ {Cy, Vp G N, Vfc G N. 

The Floquet multipliers of the multi-periodic recurrence fITS]) . are the two 
roots of the quadratic equation 

\^-{Tr C)A + detC = 0. 


It is easy to see that det C 


f{T-l) 

/(-I) 


T-l 


i=o 
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